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The fermion sign problem appearing in the mean-field approximation is considered, and the sys¬ 
tematic computational scheme of the free energy is devised by using the Lefschetz-thimble method. 
We show that the Lefschetz-thimble method respects the reflection symmetry, which makes physical 
quantities manifestly real at any order of approximations using complex saddle points. The formula 
is demonstrated through the Airy integral as an example, and its application to the Polyakov-loop 
effective model of dense QCD is discussed in detail. 
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I. INTRODUCTION 

The fermion sign problem is a big obstacle for non- 
perturbative studies of various quantum systems rang¬ 
ing from condensed matter [1] to hadron and particle 
physics [2]. The Boltzmann weight in the path-integral 
expression becomes oscillatory, and importance sampling 
of Monte Carlo simulations fails [3] . 

Finite-temperature quantum chromodynamics (QCD) 
has been well understood due to the success of lattice 
Monte Carlo simulations [4] and of effective-model de¬ 
scriptions of the confinement-deconfinement crossover [5- 
10]. On the other hand, a lattice simulation of dense 
QCD faces the fermion sign problem [2], which prevents 
us from exploring properties of QCD matters under ex¬ 
treme conditions, e.g., in colliders and in neutron stars 
[11, 12]. Moreover, it causes a severe disease in model 
studies as well: The effective action as a function of order 
parameters takes complex values [13, 14]. The evaluation 
of the mean-field free energy requires an oscillatory inte¬ 
gral, and the sign problem remains [15]. This situation 
motivates us to devise a systematic and exact computa¬ 
tional scheme that evades the sign problem. 

Recently, the Picard-Lefschetz theory turns out to be 
helpful for evaluating such oscillatory integrations by us¬ 
ing the saddle-point analysis in the complexified space of 
field variables [16-21], and it is applied to lattice Monte 
Carlo simulation [22-26]. Other progress has been made 
in the Polyakov-loop effective model: Even after the ana¬ 
lytic continuation of gauge fields, the free energy becomes 
real by choosing a complex saddle point that respects 
charge- and complex-conjugation symmetry [27, 28]. 

Unifying these ideas of recent progress, we derive a 
mathematical formula that allows us to systematically 
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calculate the partition function. The formula is especially 
suitable for the saddle-point analysis, and we show that 
the formula makes the partition function real valued at 
any order of approximations using complex saddle points. 
We demonstrate its usefulness through an example of the 
Airy integral. Its application to QCD is also discussed, 
and we show the result of the saddle-point approximation 
in the SU (3) matrix model of heavy dense quarks. 


II. FORMALISM AND MAIN THEOREM 

Let us consider a multiple integration that gives the 
partition function, 

Z= f (1) 

where S{x) is the action functional of the real field x = 
(xi,... ,Xn)- By definition, the partition function (1) 
for physical systems must be a real quantity, however 
the Boltzmann weight S{x) can be complex in general. 
One of the sufficient conditions for ensuring Z £ K is the 
existence of a reflection symmetry L : (xi) i-^ (LijXj), 
which satisfies Lij = Lji £ K, = 1 and 

W^=S{L-x). ( 2 ) 

In the case of finite-density QCD, we will see later that 
the charge conjugation gives this correspondence. The 
main purpose of this section is to derive the useful for¬ 
mula for the mean-field analysis and its systematic im¬ 
provement while keeping Z real under Eq.(2). 

For evaluating oscillatory multiple integrations, the 
Picard-Lefschetz theory is powerful and its application 
to the path integral has recently gathered much atten¬ 
tion [16-26]: This method complexifies the space of field 
variables and finds the most convergent integration cy¬ 
cle around the saddle point, which is called the Lefschetz 
thimble. Let us denote complexified variables of a; as z. 
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To construct the Lefschetz thimble, we look for the 
saddle points z'^ in the complexified field space, which 
are formally labeled by cr S S. The directions of quan¬ 
tum fluctuation around 0 °’ can be designated by solving 
Morse’s flow equation [16-19]: 


dz^ _ / dS{z) \ 
dt V 9zi ) 


( 3 ) 


The Lefschetz thimble Za is identified as the set of points 
reached by some flows emanating from z'^. As a conse¬ 
quence, each Z(T forms an n-dimensional real manifold. 
The partition function can now be computed as 

Z='^n^ [ d”z (4) 

(tGS 

Therefore, the original integration (1) is now decomposed 
into the sum of the nicely converging integrations around 
the saddle points z'^. The coefficient is given by the 
(oriented) intersection number between K" and = 

(Acr,R"). The dual thimble is defined as the set of 
the points reached by flows getting sucked into 2 °’. 

Let us rewrite Eq.(4) into the formula manifesting Z € 
K. We extend the linear map L to an antilinear map 


K : (zi) ^ (LijZj). (5) 

The flow equation (3) shows the covariance under this 
antilinear reflection: Using Eq.(2), we find 

<«) 

and thus the antilinearly transformed function K{z(t)) 
satisfies the same flow equation (3). Many important 
properties can be deduced from this conclusion: 


respondence between S_|_ and E_. Equation (4) becomes 



Each integral on the r.h.s. of the formula (8) is real 
or purely imaginary depending on whether K changes 
orientation of Jo- U ■ Since the l.h.s. is real, rir must 
be zero unless the integral on J,- -|- Zr is real [30] . This 
conclusion can also be applied to expectation values of 
any physical observables that satisfy the symmetry (2). 
The decomposition formula (8) takes a suitable form for 
the saddle-point analysis. 

Let us emphasize here that the reality of the partition 
function is ensured only by the invariance under the an¬ 
tilinear map K. The point of our discussion is that the 
Lefschetz-thimble decomposition of the integration cycle 
manifestly respects the antilinear reflection K, and so 
does the saddle-point analysis based on it. 

There are two possible remnants of the sign problem 
in the Lefschetz-thimble method [24, 25, 31, 32]. One is 
the complex phases coming from the Jacobian of d^z, 
which is called the residual sign problem. In Eq.(8), 
the imaginary part of the Jacobian automatically can¬ 
cels under the antilinear symmetry AT, but its real part 
can flip the sign and gives the sign problem when one 
goes beyond the saddle-point approximation. Another 
one comes from the overall signature of each integral in 
Eq.(8), or the signature of rig.. A cancellation among 
different saddle-point contributions is important for the 
study of the phase transition [21], because the partition 
function vanishes at the phase transition [33]. This can 
give information about the phase transition even within 
the saddle-point analysis. 


1. When z’^ is a saddle point, so is K{z^). Let 

be the Lefschetz thimble around z'’"; then Za '■= 
{K{z) I z e Jo-} coincides with the Lefschetz thim¬ 
ble around K{z'^) up to orientation. 

2. Since the flow is covariant under the antilinear map, 

(Acr,R") = under an appropriate choice 

of the orientation. This means that the contribut¬ 
ing thimbles always form an invariant pair Za^Za ■ 

Let us decompose E into three disjoint parts. For sim¬ 
plicity, we assume that S{z'^) S K only if the saddle point 
satisfies z°’ = K{z'^) [29]; then, E = EoUE+UE_, where 

Eo = {cr I z'^ = L • P"}, 

E+ = {a I Im5'(z'^) > 0}, 

E_ = {ct I Im5'(z'") < 0}. (7) 

The transformation K induces a one-to-one and onto cor- 


III. EXAMPLE: AIRY INTEGRAL 

Let us work on the simplest example in order to un¬ 
derstand the formula (8). We consider the Airy integral 
with a real parameter a: 

zAiry {a) = [ dx (9) 

Denoting the exponent of the integrand as —S{x,a), 
S{x,a) = S{—x,a) when x and a are real, and (2) is 
satisfied. The antilinear map K is then K{z) = —z. 

If a > 0, then the saddle points are z"^ = ±-\/ai, and 
both of them are invariant under K; z'^ = —z”. Study¬ 
ing the flow equation, only one saddle point +\/ai turns 
out to contribute to Eq.(9). Using the Lefschetz-thimble 
decomposition (8), the saddle-point approximation gives 

^Airy(a) ~ exp - (10) 
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which is the asymptotic formula of the Airy function as 
a —>■ +00 (see Chap. 19 of [34]). On the other hand, when 
a < 0, the saddle points form the conjugate pair under 
K, ±Y^|a[, and both of them give contribution to the 
original integration. Similarly, we get 


sentation 3 is an order parameter for the confinement- 
deconfinement transition [35]; 


4 = 



V exp 



(16) 


■^Airy (a) 





( 11 ) 


which is nothing but the asymptotic formula of the Airy 
function as a —>■ —oo (see Chap.19 of [34]). Notice that 
both approximate results [Eqs. (10) and (11)] are mani¬ 
festly real, as we have already shown in general. 

Let us remark that only by using S{x,a) = S{—x,a) 

is the partition function ensured to be real: 

1 fx^ \ 

■^Airy = dx COS - \- GX ] . (12) 

Jo V 3 J 


where V refers to the path ordering. To understand its 
properties, the constrained effective action 5 'eff(a 4 ) of the 
Polyakov loop is useful [36, 37]. It describes the (com¬ 
plex) weight of the partition function under the back¬ 
ground Polyakov-loop phase 04 : 

exp(-S'eff(a 4 )) = y PA ( 5(^4 - 04 ), (17) 

where A 4 is the temporal gauge field and S'[A] = 
Sym[A] — lndetAI(/rqk, A). We implicitly take the 
Polyakov gauge fixing with 


However, this integration is oscillatory and still has the 
sign problem. The Lefschetz-thimble decomposition (8) 
plays a pivotal role in evaluating these asymptotic formu¬ 
las (10) and (11) without encountering the sign problem. 


IV. APPLICATION TO THE SIGN 
PROBLEM OF DENSE QCD 


The QCD partition function at temperature T = jS ^ 
and quark chemical potential /Tqk is given by 

^QCD = / detM(^qk, A) (13) 


where Sym = ytr Jq da:"^ / d^x\F^^\‘^ (> 0) is the Yang- 
Mills action, and 

detAI(/rqk, A) = det + i^A^) + y'^/Tqk + ruqk] 

(14) 

is the quark determinant. When / 4 qk 7 ^ 0, the quark de¬ 
terminant becomes an oscillatory functional of the gauge 
field A, and the sign problem emerges. Even when 
^qk 7 ^ 0, the charge conjugation A 1-7 —A* with the 75 
hermiticity implies that the fermion determinant still sat¬ 
isfies the identity [14, 15], 

det AI(/rqk, A) = det AI(—/Tqk, A'^) 

= det AI(/rqk,-A). (15) 

The charge C and complex 1C conjugation, or the C/C 
transformation, serves as the antilinear map (5) for finite- 
density QCD [27, 28]. The Lefschetz-thimble decompo¬ 
sition ( 8 ) manifestly respects the C/C symmetry so that 
^QCD G K. Since the discussion is robust, this conclusion 
applies to any effective descriptions of QCD, including 
lattice QCD simulations [22-26] and also matrix models. 

We apply our insight on the Lefschetz-thimble method 
to the sign problem in Polyakov-loop effective mod¬ 
els. The Polyakov loop £3 of the fundamental repre¬ 


exp(i^) = diag 


gi(0i+02) gi( —^1+02) g—2i02 


(18) 


The Weyl group action ( 6 * 1 , 02 ) (— 0 i, 02 ), ( 0 i, 02 ) 

((01 + 362 )/ 2 , (01 — 02 )/ 2 ) just permutes eigenvalues of 
the Polyakov loop (18), and thus the parameter region 
can be restricted to £ = {0 = ( 0 i, 02 ) | 3|02| < 0i < tt}. 
This background field method, or the mean-field approx¬ 
imation (MFA), gives a single integration over 04 ( 01 , 02 ) 
to compute the partition function: 


Zqcd — [ d0id02i7(0i, 02) exp [— S'eff(0i, 02)], (19) 

with H{9) = sin^ 0i sin^((0i + 302)/2) sin^((0i — 302)/2) 
coming from the SU{3) Haar measure [7, 9]. When 
the quark chemical potential /Xqk is turned on under the 
nontrivial Polyakov-loop background, the effective action 
Ses{0) takes complex values in general due to the quark 
determinant. This makes the integration (19) oscillatory, 
and the sign problem remains in the MFA [14, 15]. 

In order to evade the sign problem of the MFA, we 
should rewrite Eq.(19) based on the decomposition for¬ 
mula (8) . Since H (0) identically vanishes at the edge of £, 
the Lefschetz-thimble method can be applied by regard¬ 
ing S{9) = 5'eff(0) — lniL(0). In terms of the Polyakov- 
loop (18), the C/C transformation K becomes 


K{zi,Z 2 ) = {zi,-Z 2 ), (20) 

where Zi is the complexified variable of 0^, up to some 
Weyl transformations. The C/C symmetry (15) leads 


Ses{zi,Z2) = Ses{zi, - 22 )- ( 21 ) 

It is important to remark that the Polyakov loop, 

4(0) = ^ (2e‘®= cos 01 + , (22) 

satisfies the C/C symmetry (21). For our demonstration. 
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FIG. 1. Flows at /i = 0.1 and fi = 2 around the C/C-symmetric 
saddle point (the black blob) z* in the Re( 2 i)-Ini( 2 ; 2 ) plane. 
The red solid and green dashed lines are the Lefschetz thimble 
3* and its dual .R* of z*, respectively. Since At, intersects with 
Im( 22 ) = 0, the integral over 3* is equal to that over £. 


1 in this limit, as was numerically observed in Fig.2 of 
Ref. [14]. Using the saddle-point approximation, we can 
see that these Polyakov loops have different expectation 
values: 

(%) ~ (^ 3 ) — T (smh2iz2 — 2cosz* sinhi^^) > 0.(25) 

O 

These results are consistent with those of exact compu¬ 
tations of expectation values [14], and the C/C-symmetric 
saddle point gives real expectation values [27, 28]. The 
formula (8) shows that this approximation can be system¬ 
atically improved perturbatively to satisfy these physical 
requirements. 


V. CONCLUSION AND DISCUSSION 


we take a simplified heavy-quark model [14, 38-43]: 

~ (e^4(^l,g2) +e-^%(gi,g2)) 

= — ^ ^2 cos 01 cos(02 — im) + cos(202 + ■ (23) 


When h ^ 0 and ^ ^ 0, the integration (19) is 

oscillatory because takes complex values. The gluon 
dynamics is neglected just for simplicity. 

If h = 0, the system has the Z 3 symmetry generated by 
( 01 , 02 ) ( 01,02 + 27r/3) up to Weyl group actions. The 

Haar measure factor 7 J( 0 i, 02 ) takes the unique maxima 
at 0* = (27r/3,0) in C, and this is Z 3 invariant. For 
later purposes, it is important to notice that this saddle 
point is CK, symmetric in the sense that 0* = K{9*). 
The eigenvalues of the temporal Wilson line (18) become 
uniformly separated, and thus (£3) = 0 . 

When 0. 7 ^ 0, the Z 3 symmetry is explicitly broken, and 
the saddle point z* moves away from the Z 3 -symmetric 
point (27r/3,0), but it still must be CK, symmetric be¬ 
cause there are no other saddles to form a pair. By conti¬ 
nuity of intersection numbers, z* will contribute to Zqcu 
even for finite h and /r. In order to demonstrate it lucidly, 
we explicitly solve the flow equation (3) around the C/C- 
symmetric saddle point (see Fig. 1). The dual thimble M.* 
of z* is shown with the green dashed curve, and it indeed 
intersects with the original integration cycle £. There¬ 
fore, the complex saddle point contributes, but physical 
quantities take real values since the CK symmetry is re¬ 
spected under the Lefschetz-thimble decomposition. 

In the limit fj, —7 -l-oo, the saddle-point approximation 
becomes accurate. In this limit, the effect of the quark 
potential S^s becomes quite large. By solving the saddle- 
point equation of S in the limit /i —>■ 00 , we find that 


z 


* 

1 


^ ,e->^ 


(24) 


and it indeed approaches the perturbative vacuum 
Zpert. = (0,0) of ^eff. Thus, (£3) and (£3) converge to 


In this paper, we have shown that the Lefschetz- 
thimble decomposition manifestly respects the CK sym¬ 
metry of dense QCD, so that Zqcd becomes real even 
with approximations. The Lefschetz-thimble decomposi¬ 
tion is suitable for the saddle-point analysis, so the parti¬ 
tion function can be computed in a systematic way. This 
solves the sign problem appearing in the MFA when the 
classical action is complex. The importance of such an- 
tilinear reflection symmetry has also been discussed in 
the literature of RT-symmetric quantum theory and the 
reality of the partition function was also ensured there 
[44-46]. It will be interesting to establish a clear and 
explicit connection between these two consistent results. 

Our method was demonstrated by applying it to the 
Airy integral and to the SU(3) matrix model of heavy 
dense quarks at finite /r. Since the derivation is robust, 
our result can be broadly applied: In condensed mat¬ 
ter physics, studying Hubbard models with Hubbard- 
Stratonovich transformation must be an interesting ex¬ 
ample [47], and in hadron physics other effective models 
of QCD—such as Polyakov-loop extended Nambu-Jona- 
Lasinio model [10] —should also be studied based on the 
formula ( 8 ). There has been empirical proposals to cir¬ 
cumvent the sign problem in the MFA [14, 15, 48, 49], 
and the reality of the free energy at the C/C-symmetric 
saddle point has been argued in Refs. [27, 28]. The cur¬ 
rent work goes beyond the saddle-point approximation 
and gives an exact expression for the partition function 
using the Lefschetz-thimble technique ( 8 ). This approach 
provides a systematic way to find relevant saddle points, 
and we now have a solid mathematical foundation for 
those approximations. This work will give clues for in¬ 
vestigating the sign problem of lattice field theories. 

We should note that, however, the positivity of the 
partition function is not yet ensured in general cases, es¬ 
pecially when many saddle points contribute to the parti¬ 
tion function. This possible cancellation among contribu¬ 
tions of several saddle points will be physically important 
in the context of phase transitions. In order to deepen 
our understanding of positivity and also of unitarity, we 
need to combine our result of this paper with other prop- 
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erties of physical systems, and we expect that it might be 
related to the reflection positivity [50] . This would be an 
important direction for future studies aimed at solving 
the sign problem of other complicated systems like chiral 
random matrix models [51]. 
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